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Abstract. We have obtained the full 1-200 spectrum 
of the low luminosity (36 L©) Class I protostar Elias 29 
in the p Ophiuchi molecular cloud. It provides a unique 
opportunity to study the origin and evolution of interstel- 
lar ice and the interrelationship of interstellar ice and hot 
core gases around low mass protostars. We see abundant 
hot CO and H2O well as the absorption bands of 

CO, CO2, H2O and "6.85 /im" ices. We compare the abun- 
dances and physical conditions of the gas and ices toward 
Elias 29 with the conditions around several well studied lu- 
minous, high mass protostars. The high gas temperature 
and gas/solid ratios resemble those of relatively evolved 
high mass objects (e.g. GL 2591). However, none of the ice 
band profiles shows evidence for significant thermal pro- 
cessing, and in this respect Elias 29 resembles the least 
evolved luminous protostars, such as NGC 7538 : IRS9. 
Thus we conclude that the heating of the envelope of the 
low mass object Elias 29 is qualitatively different from 
that of high mass protostars. This is possibly related to a 
different density gradient of the envelope or shielding of 
the ices in a circumstellar disk. This result is important 
for our understanding of the evolution of interstellar ices, 
and their relation to cometary ices. 

Key words: Stars: formation - individual: Elias 29 - ISM: 
dust, extinction - molecules - abundances - Infrared: ISM: 
lines and bands 
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1. Introduction 



The general picture of low mass star formation has 
been formed since the 1980's with the availability of in- 
frared and millimeter wavelength broad band photometry 
from the ground, and with the IR AS sa tellite and KAOob- 
servatory (e.g., La da fc Wilking 1984; Adams et al. 1987; 

1993| ). A classifica- 



Hillenbrand et al. 1992; Andre et al 



tion scheme was made, where the continuum emission of 
Class and I objects peaks in the submillimeter and far- 
infrared. These objects are still deeply embedded in their 
accreting envelopes. In the Class II phase, the wind of the 
protostar has cleared its surrounding environment, such 
that it becomes optically visible, and shows Hi emission 
lines. The continuum emission of these objects peaks in 
the near-infrared, but there is still significant excess emis- 
sion above the stellar continuum. They are believed to be 
surrounded by optically thick dusty disks. Finally, little 
dust emission remains for Class HI objects, when the disk 
is optically thin, and planetary companions may have been 
formed. 

Our knowledge of the physical and chemical state and 
evolution of the material surrounding protostars, has pro- 
gressed with the availability of medium and high res- 
olution spectroscopic instrumentation at near and mid- 
infrared wavelengths (^^ 2 — 20 /im). The progress made, 
is best illustrated by the observations of high mass 
protostars, which are bright and easy to observe. The 
(ro-)vibrational bands of various molecules (CO, H2O, 
CH3OH, silicates) were observed from the ground, re- 
vealing profile variations of the 3.07 /im, and 4.67 /im 
H2O and C O ice bands (e.g. Smith et al. |1989| ; Tielens 

199^ ). This was interpreted as 



et al. 1991; Chiar et al. 



evaporation of the volatile CO ice, and crystallization of 
H2O ice in the molecular envelopes. This heating effect is 
strengthened by the detection of hot gas in 4.6 /im CO 
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observations (Mitchell et al. |l99l| ). With the launch of 
the Infra red Space Observatory in 1995 (ISO; Kessler et 
al. 1996), it became possible to observe all other molec- 
ular bands in the infrared (Whittet et al. 19961 ). It was 
shown that high mass protostellar evolution can be traced 
in the gas-to-solid abundance ratios (van Dishoeck et al. 



1996; van Dishoeck & Blake 1998), and the profiles of the 



199S 



ice bands, in particular solid CO2 (Gerakines et al 
Boogert et al. 2000| ). Thus, there is overwhelming evidence 
that thermal processing, i.e. evaporation and crystalliza- 
tion of ices in and around hot molecular cores, plays an 
important role in the evolution of high mass molecular 
envelopes. 

The composition and evolution of the molecular ma- 
terial around low mass protostars are not as well studied. 
It seems unlikely that the molecular material evolves sim- 
ilar to that around high mass protostars. Low mass pro- 
tostars evolve much slower, release less radiative energy, 
drive less energetic winds, and form disks. It is not estab- 
lished whether low mass objects possess hot cores as well, 
and whether the ices survive the process of star formation. 
If (some of) the ices survive, are they included into comets, 
and if so, are the ice structure and composition still the 
same compared to interstellar ices? How important are 
energetic processes, such as cosmic ray bombardment, in 
altering the ice composition on the long time scale of the 
formation of low mass stars? 

To investigate the influence of low mass protostars on 
their molecular envelope, we make an infrared spectro- 
scopic study Elias 29, also called WL 15 and YLW 7 (Elias 

1986| ). On a large 



197S, Wilking et al. 1983, Young et al 



scale, Elias 29 lies in core E, which is in the south-east 
corner of the 1x2 pc extended compact CO ridge L 1688 



(Loren et al. 1990) in the densest part of the p Ophiuchi 
cloud, at a distance of ~160 pc from the earth (Wilking & 



Lada[1983[ Whittet |1974| ). It is the reddest object foun d 
in the near- infrared survey of this cloud by Elias ( 1978|) , 
without a counterpart at optical wavelengths. For our ob- 
servations, we used Elias' coordinates (J2000): 



a = 16^27™09".3 



S = -24°37'21' 



The overall spectrum of Elias 29 is typical for a heav- 
ily embedded Class I s ource , probably in a late accretion 
phase (Wilkin g et a l. |1989| ; Andre fc M ontmerle |l994 ; 
Greene & Lada 1996| ; Saraceno et al. 199(j ). The embedded 
nature is also revealed by its high extinction, and by the 
cold compact envelope observed at millimeter wavelengths 
(Andre & Montmerle |1994| Motte et al. |199S| ). Ehas 29 
is associated with a molecular outflow (Bontemps et al. 
1996| ; Sekimoto et al. 1997). With a bolometric luminos- 
ity of ~ 36 Lq (Chen et al. |1995D, Ehas 29 is the most 



luminous protostar in the p Oph cloud, which makes this 
source very suitable for spectroscopic studies. The rela- 
tively high luminosity, and high bolometric temperature 



(Chen et al. 1995). In the pr e-mai n-sequence evolution- 
ary tracks of Palla & Stabler (1993), this corresponds to 
a star with end mass 3.0-3.5 Mq. Elias 29 might thus be 
a precursor Herbig AeBe star. This classification is how- 
ever uncertain. For example, it has been argued from the 
SED, and the absence of mid-infrared emission features, 
that Elias 29 is a 1 Mq protostar with a large accretion 
luminosity, and a spectral type of K3-4 at the birth line 



(Greene & Lada ^OOOj ). 

This Paper is structured as follows. Technical details 
on the ISO infrared observations are given in Sect. ||. All 
the observed emission and absorption features are dis- 
cussed in detail in Sect. ||. Section 3J gives a description 
of the continuum shape, and a comparison to other lines 
of sight. The ice composition and thermal history, and the 
silicate band depth with inferred extinction and column 



densities toward Elias 29 are discussed in Sect. 3.2. Then 



numerous lines of gaseous CO and II2O are detected, and 
modeled to derive gas temperatures and column densities 



(Sect. 3.3). The molecular abundances and gas-to-solid ra- 
tios of Elias 29 are compared to a sample of sight-lines, 
ranging from dark cloud cores to evolved protostars. A 



compar ison with high mass protostars is made (Sect. 4.1). 
Section L2 discusses the origin of the wealth of observed 
emission and absorption features and puts them in a ge- 
ometrical picture, where we review the evidence for an 
extended envelope and an accretion disk. We conclude in 
Sect. ^ with a summary and suggestions for future obser- 
vations. 



2. Observations 



2.1. The 2.3-45 p,m spectrum 



A low resolution {R = A/AA = 400), fuU 2.3-45 iim 
spectrum of Elias 29 was obtained with the ISO Short 
Wavelength Spectrometer (ISO-SWS; de Graauw et al. 



19961 ) during revolution 267 (August 10 1996). The ISO- 



bol 



410 K) imply an age in the range 0.5-4 10 yr 



SWS pipeline and calibration files, available in July 1998 
at SRON Groningen, were applied. The spectrum is gen- 
erally of good quality, with well-matching up and down 
scans, and no serious dark current problems, except for 
band 2C (7-12 /im). Here, we found that the up and down 
scans deviate over the silicate band. One scan showed 
good agreement with a ground-based spectrum of Han- 
ner et al. ( |1995D , and we used this to correct the deviating 
scan. Standard after-pipeline steps were applied, such as 
low order flat-fielding, sigma clipping and re-binning (see 
also Boogert et al. 199S| ). The twelve sub-spectra in the 
2-45 /im range match fairly well at the overlap regions. 
Small correction factors (<15%) were applied to correct 
for the band jumps. 



3 



At selected wavelength ranges (3-3.6, 4-9, and 19.5- 
28 A*m), we also obtained high resolution {R = 1500) 
ISO-SWS grating spectra, in revolution 292 (September 
04 1996). These were reduced similarly to the low reso- 
lution spectrum. We found that the overall shape of the 
spectrum near 4-5 /im is quite badly affected by detector 
memory effects, presumably due to the occurrence of scan 
breaks (de Graauw et al. 1996| ). We corrected for this, 
by applying a wavelength-dependent shift to match the 
low resolution spectrum. This does not affect our conclu- 
sions, since the high resolution spectrum was only used to 
study narrow features. Also, near 6.9 fj,m the scans devi- 
ate significantly because of memory effects. This problem 
is reflected in the large error bars given in this paper, as 
they were derived from the difference between the average 
up and down scans. 

2.2. The 45-190 ijm spectrum 

Elias 29 was observed during Revolution 484 (March 14 
1997) with the ISO Lon g Wavelength Spectrometer (ISO- 
LWS; Clegg et al. 1996 ). We obtained 15 scans covering 
the range from 43 fim to 197 /zm in the low resolution 
mode {R ^--^200) for a total of 2611 sec of integration 
time. The data were reduced using the Off-Line-Processing 
package (OLP) version 7 and the ISO-Spectral-Analysis- 
Package (ISAP) version 1.3. The sp ectra were flux cali- 
brated using Uranus (Swinyard et al. 1996 ). We find that 
at the ISO-SWS/LWS overlap region near 45 /zm, LWS 
has 35% higher flux than SWS. This difference is only 
slightly larger than the absolute calibration uncertainties 
of the two instruments, and thus it is doubtful that this 
can be ascribed to the presence of extended emission in 
the larger aperture of ISO-LWS (~ 80" versus ~ 25"). 
We therefore decided to multiply the LWS spectrum down 
with this factor. 

3. Results 



3.1. The spectral energy distribution (SED) 

Elias 29 is only visible at wavelengths larger than 
~1.5 /xm (Greene & Lada |1996| ; Elias |1978D . Our ISO ob- 
servations show that the continuum emission rises steeply 
between 2-3 /xm, reaches a maximum of Ai^A=15xlO~^^ 
W cm~^ at A ^--^5 /im, and is remarkably flat with \Fx^ 
8 X 10-16 W cm-2 between 20 and 100 ^lui (Fig. g). The 
emission has dropped to \F\^ 4 x 10-^^ W cm"^ at 
200 /xm, and by four orders of magnitude at 1300 /im. Our 
near-infrared spectral continuum fluxes are in excellent 
agreement with broad band fluxes from ground-based ob- 
servations (Elias 197^ ) . Also the ground-based small beam 



10 and 20 /im observations, as well as the large beam 12 
and 25 /im IRAS fluxes, match the ISO-SWS observa- 
tion well, thus indicating that at these wavelengths the 
emission is well confined within a region of 8" in diam- 
eter (Fig. 0; Lada & Wilking |1984| ; Young et al. |1986D . 



The reasonable match of the ISO-SWS and LWS spectra 
(Sect. 2.2) indicates that also at 45 /im the emission is not 
very extended (< 25"). At 100 /im, however, some large 
scale emission may be present, since the IRAS flux, ob- 
served in a 5.5 times larger aperture, is a factor of 2 larger 
compared to the ISO measurement. 

The observed SED of Elias 29 is different from that 
of massive protostars such as GL 2591, and GL 7009S, 
which peak in the far-infrared (Fig. It has been pro- 
posed that the shape of SEDs is independent of the lu- 
minosity of the central object, and is rather determined 
by the total dust column density along the line of sight 
(Ivezic & Elitzur 1997). In the "standard model" of Ivezic 
& Ehtzur, GL 2591 and GL 7009S would have a column 
density corresponding to an Ay of several hundred. Elias 
29 must have an Ay < 100, because it does not peak in 
the far infrared. However, the flatness of the SED up to 
100 fiia is not reproduced in these models. A lower column 
density alone thus cannot explain the differences between 
the SED of Elias 29 and massive protostars. Other factors, 
such as a different density gradient, and the presence of a 
circumstellar disk are probably important. We will discuss 
the structure of Elias 29, in relation to the detected gas 



and ice absorption features, in Sect. 4.2 



3.2. Ice and dust absorption bands 

Numerous absorption bands of ices and silicates are 
present in the infrared spectrum of Elias 29 (Fig. ||), such 
that hardly any 'clean' continuum emission is left. We 
identify each band, derive column densities, and, when 
possible, determine the ice mantle composition and ther- 
mal history. The full spectrum also allows a determination 
of upper limits of abundances for undetected, though as- 
trophysically relevant molecules. The column densities for 
the species discussed below are summarized in Table |l|. 

3.2.1. H2O ice 

The infrared spectrum of Elias 29 shows all the vibration 
modes of H2O ice in absorption (Fig. ^). We see the O- 
H stretching mode at 3.0 /im ('Vi, 1^3" in spectroscopic 
notation), the 0-H bending mode at 6.0 /im ("^2"), the 
libration or hindered rotation mode at ~12 /im ('Vl"), 
the combination mode at 4.5 /im ("3^'l" or 'V2 -I- J^l"); 
and perhaps the lattice mode at ~45 /im. 

The continuum determination is complicated by the 
large width of all these bands. In accordance with other 
studies (Smith et al. 1989| ; Schutte et al. 1996; Kcane et 



al. 2000), we used single blackbodies to fit the continuum 
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Fig. 1. Spectral energy distribution 
of Elias 29, consisting of ground- 
based observations (A <2.4 /im; 
Greene & Lada |l99^), an ISO- 
SWS spectrum (A =2.4-45 ^ni), 
and an ISO-LWS spectrum (A =45- 
195 /im). The data point at 1300 /im 
i s tak en from Andre & Montmerle 
( 1994 ), which we have connected 
with a dashed straight line to the 
1S0~LWS spectrum, to guide the eye. 
The dotted line is the adopted con- 
tinuum, as determined by blackbody 
fits and by hand. The open circles 
are ground-based and IRAS observa- 
tions (see text). The top inset shows 
a magnification of the 10-200 /im re- 
gion. 
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Fig. 2. Spectral energy distribution 
of Elias 29 compared to the high 
mass protos tars GL 7009S (Dartois 
et al. |l998a|), a nd GL 2591 (van der 
Tak et al. |l999| ). The spectrum of the 
Herbig Ac star AB Aur is a compila- 
tion of continuum observations taken 



from Mannings (1994), and is further 
discussed in Sect. 1.2. The flux scale 



of each spectrum has been divided by 
the values given in brackets. 



locally, directly adjacent to the absorption bands. For the 
6 /im band we took into account that laboratory spectra 
of the bending mode of H2O ice show a prominent wing on 
the long wavelength side, extended up to 8 /im (e.g. Hud- 
gins et al. 1993, Maldoni et al. 1998). We simultaneously 



fitted a blackbody continuum, normalized at 5.1 /im, and 



a laboratory ice spectrum to the observed flux at 8 /im 
and the shape of the 6.0 /im feature (Fig. |). 

The shape of the 6.0 /im H2O bending mode is partic- 
ularly sensitive to the ice temperature (e.g. Maldoni et al. 
1998). At higher T in the laboratory, the strength of the 
main 6.0 /im component decreases at the expense of more 
absorption in the long wavelength wing. In the spectrum 
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Fig. 3. Low resolution (i? = 400) mid-infrared ISO-SWS spectrum of Elias 29 with a smooth, global continuum 
(dotted line), rather arbitrarily determined by hand and blackbody fits. The vibrational absorption bands of various 
molecules are indicated. 



of Elias 29, the wing cannot be seen as a separate feature, 
since at 8 /im it blends with the very deep silicate band. 
However, the observed 6.0 /im band is relatively sharp, 
and it can only be fitted with H2O ice at T < 80 K, with 
a best fit at r = 40 K (Fig. ||). The excellent fit to the 
6.0 /im band in Elias 29 indicates that the 5.83 and 6.2 /im 
excess absorptions detected toward several massive proto- 
stars (Schutte et al. |1996| , |1998| ; Keane et al. |2000D , are 
not seen in this source (Sect. 3.2.7). 

The observed peak position of the stretching mode 
of H2O ice toward Elias 29 is 3.07±0.01 [xm. The short 
wavelength wing is well matched with a laboratory ice at 
T =40 K, as for the bending mode (Fig. ^). The long wave- 
length wing, however, is poorly fitted. It has been realized 
since long that light scattering by large ice grains leads to 



Table 1. Column densities of ices toward Elias 29. Non- 
detections are indicated with 3cr upper limits. 



molecule 


N [10^^ 


H2O 


34 (6) 


"CO2 


6.7 (0.5) 


"CO2 


0.083 (0.005) 


CO 


1.7 (0.3) 


CH4 


< 0.5 


NH3 


< 3.5 


CH3OH 


< 1.5 


H2CO 


< 0.6 


HCOOH 


< 0.3 


OCS 


< 0.015 


XCN 


< 0.067 



extra extinction on the long wavelength wing (e.g. Leger 
et al. |1983D . To illustrate this, we calculate the extinc- 
tion cross section for spherical silicate grains coated with 
ice mantles, applying the code given in Bohren & Huff- 



man (1983) and the optical constants of Draine & Lee 
(1984) and Hudgins et al. (1993). Indeed, grains with a 
core+mantle radius of ^0.6 /im provide a much better fit 
to the long wavelength wing than small grains do (Fig. H). 
This effect is unimportant for the 6.0 /im band since it is 
intrinsically weaker, and the grains are smaller compared 
to the wavelength. In a more realistic approach, a distribu- 
tion of grain sizes, as well as constraints to other observ- 
ables such as continuum extinction, the total grain and 
ice column densities, and polarization need to be taken 
into account. Although there is a general consensu s that 
large grains need to be invoked (e.g. Leger et al. 1983; 



Pe ndleto n et al. |1990| ; Smith et al. |1989| ; Martin & Whit- 
tet 1990 ), there is no unified grain model yet tha t obe ys all 
the observational constraints (e.g. Smith et al. 1993; Tie- 
lens 1982). Alternative absorbers at the long wavelength 



wing have been proposed, such as H2O.NH3 bondings. As 
illustrated in Fig. this effect is however small at the low 
column density ratio of NH3/H2O < 0.13 toward Elias 29 
(Sect. 3.2.7). A small contribution is also made by absorp- 
tion by hydrocarbons (Sect. 3.2.4). 

The peak optical depth of the 3.0 /im band is 
1.85±0.08, which is in excellent agreement with the study 



of Tanaka et al. (1990). Using an integrated band strength 
A = 2.0 X 10^^® cm molecule"^, we derive a column den- 
sity of iV(H20)= 3.0 X lO^** cm-2 for the small grain 
model, and 3.7 x 10^* cm^^ for the large grain model. 
Since at present we can not favor one of these two cases, we 
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Fig. 4. Optical depth spectrum of Elias 29, assuming the 
continuum indicated in Figs. |l] and ^ The light, thick line 
is a labor atory spectrum of H2O ice at T =10 K (Hud- 
gins et al. 1993). All five H2O ice vibration bands can be 



discerned. 
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will assume an average value of iV(H20)= (3.4±0.6) x 10 
cm~^ in this paper. The error bar also includes the uncer- 
tainty in band strength, which increases with 10% when 
the ice is heated from 10 to 100 K (Gerakines et al. 1995| ). 
Note that a column density determination from the 6.0 /im 
bending mode is more uncertain due to the unreliable con- 
tinuum on the long wavelength side (Fig. ||). At this col- 
umn density of H2O ice, the depth of the other vibrational 
modes is in good agreement with the observed spectrum 
of Elias 29 (Fig. |). 



3.2.2. CO ice 

The CO ice band at 4.67 /im in Elias 29 is contami- 
nated by gas phase CO lines from low J levels (Figs. || 
and |l^). In particular, the P(l) line lies in the cen- 
ter of the ice band at 4.674 /im. To study the band 
profile, we subtracted a model for the gaseous lines at 
Tex = 750 K, = 5 X 10^^ cm"2, and 5d=5 km s^^ 
(Sect. ^T^ ). This increases the ice band width by 0.9 cm^^, 
to FWHM=4.40 cm-i (0.010 /im). With a peak position 
of 4.673 /im (2140.1 cm^^), the CO ice band observed to- 
ward Elias 29 is similar to that of the luminous protostar 
N GC 75 38 : IRS9 (Fig. " 
al. 



1998) 



Tielens et al. 1991; Chiar et 



The main, narrow component at 4.673 /im is at- 
tributed to pure solid CO, or CO embedded in an environ- 
ment of apolar molecules. In particular, mixtures with O2, 



at an O2/CO ratio as much as 5 (Elsila et al. 1997; Chiar et 
al. 1998) provide good fits. Mixtures of CO with CO2 are 



generally too broad (Ehrenfreund et al. 1997). While the 
apolar, volatile component dominates the spectrum, both 
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Fig. 5. Analysis of the 3.0 /im (a) and 6.0 /im (b) absorp- 
tion bands of II2O ice. In panel a, the thick gray line is 
a laboratory spectrum of pure II2O ice at T —40 K. The 
dashed line is a spectrum of H20:NH3 = 10:1 (T = 50 K). 
The dotted line gives a calculated band profile of a large 
ice grain (see text). Panel b shows the 5-8 /im spectrum, 
with laboratory spectra at T = 40 K (thick gray line), 
and at T = 100 K (thin, solid line), showing that only low 
temperatures provide good fits to Elias 29. The dashed 
line shows the assumed blackbody continuum. The nar- 
row absorption lines in the observed spectrum originate 
from II2O vapor (Sect. |3.3|). 



Elias 29 and NGC 7538 : IRS9 show evidence for a wing on 
the long wavelength side. This is attributed to CO diluted 
in a mixture of polar molecules such as H2O and CII3OH 



(Chiar et al. 1996; Tielens et al. 1991). Assuming a band 
strength A = l.lx 10^^ cm molecule ~^ for both the polar 
and apolar components (Gerakines et al. 1995| ), we derive 



N{CO ice)=1.7xl0 cm~ with an apolar/polar ratio of 
~8, comparable to NGC 7538 : IRS9. These results are in 
good agreement with the ground-based study of Kerr et 
al. ( |1993D . Although NGC 7538 : IRS9 seems to have a 
larger polar CO component in Fig. ^, this difference may 
merely reflect uncertainties in the continuum subtraction, 
and the fact that the NGC 7538 : IRS9 spectrum is not 
corrected for gas phase CO lines. 

3.2.3. CO2 ice 

The absorption bands of CO2 ice are prominently present 
in the infrared spectrum of Elias 29 (Fig. ^). We see the 
stretching and bending modes at 4.27 and 15.2 /im respec- 
tively. Not visible in this spectrum is the stretching mode 
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Fig. 6. The CO ice band on optical depth scale observed 
toward Elias 29 (thin solid line). The thick gray line is 
the spectrum with a gas phase CO model subtracted (see 
text). The dotted line represents the spectrum of the high 
mass protostar NGC 7538 : IRS9 (divided by 7.5; Tielens 



et al. 1991) showing the similarity of the band profiles. 



of solid ^^C02 at 4.38 /im, although the high resolution 
spectrum (Fig. shows a hint of its presence. A very 
sensi tive observation is presented elsewhere (Boogert et al. 
2000|) . The ^^002 bending mode and the "CO2 stretching 
mode have proven to be very sensitive to ice mantle com- 
position and thermal history. In Elias 29, these bands do 
not show the narrow substructures seen in many other pro- 
tostar s, and attributed to heat ed po lar CO2 ices (Boogert 
et al. |00C 
band (Fi, 



2000I; C 

Fiill), 



Gerakines et al. 1999). As for the CO ice 
the width and peak position of the ^^C02 
band very much resemble that of the luminous protostar 
NGC 7538 : IRS9. Thus, the CO2 ice toward Ehas 29 is 
mixed in with polar molecules, and is not much affected 
by heating. The ^^C02 column density is 22±4% relative 
to H2O ice, which is comparable to the values reported for 
high mass protostars (Gerakines et al. 1999). Finally, we 
derive an isotope ratio of ^^C02/^'^C02=81±11 in the ice 
toward Elias 29, which is well within the range found for 
the local ISM (Boogert et al. 



2000) 



3.2.4. The 3.47 ^^m band 

The long wavelength wing of the deep 3.0 /zm absorption 
band shows a change of slope at 3.38 /zm, indicative of 
a shallow absorption feature (Fig. |^). This feature is also 
detected in an independent ground based study of Elias 29 



(Brooke et al. 199£). For consistency with ground based 
studies, the continuum on each side of the feature was as- 
sumed to start at 3.37 and 3.61 //m. It must be emphasized 
however, that in particular on the long wavelength side. 




Fig. 7. a and b. Spectral structure in the long wavelength 
wing of the 3.0 /im band, a the merged high (i? — 1500) 
and low {R — 400) resolution spectra and the assumed 
polynomial continuum (dotted line), b optical depth plot 
of the detected 3.47 /im feature. The gray, thick line repre- 
sents the ground based spectrum of the high mass proto- 
starNGC 7538 : IRS9, divided by a factor of 2.3 (Brooke et 
al. |1999D , showing the C-H stretch mode of sohd CII3OH 
at 3.54 /im. This feature is absent in the spectrum of 
Ehas 29. 



the continuum is poorly defined. Fitting a smooth 6-th 
order polynomial results in an absorption band centered 
on 3.49±0.03 /xm with a peak optical depth of r=0.06 
(Fig. 0). The width is FWHM=120±40 cm-\ where the 
uncertainty includes the poorly constrained continuum on 
the long wavelength side. Features of similar width and 
peak position have been detected in several massive pro- 



tostellar objects (Allamandola et al. 1992) and in low mass 
objects and quiescent molecular cloud material (Chiar et 
al. 1996). A likely candidate for this 3.47 /im band is the 



C-H stretching mode of hydrocarbons. From the correla- 
tion of peak optical depths of this feature and the 3.0 /im 
ice band, it is concluded that the carrier for the 3.47 /im 
band resides in ices rather than in refractory dust (Brooke 
et al. |1996D . We find that with t(3.47 /im)=0.06 and 
t(3.0 /im)=1.85, Elias 29 follows this correlation very well. 

3.2.5. CH3OH ice 

In several high mass protostars, the 3.47 /im band 
is blended with a distinct narrow feature centered on 
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Fig. 8. Optical depth plot of the 5-8 fim region of Elias 29, 
after subtraction of an H2O ice spectrum at T =40 K 
as well as an H2O gas model at Tex = 300 K, iV = 
2 X 10^^ cm~^, 6d=2.5 km s~-^. This figure highlights the 
6.85 /Ltm absorption feature. The thick gray lines give a 
comparison with the massive protostars NGC 7538 : IRS9 
(top) and S 140 : IRSl (bottom; shifted down 0.08 along 
the optical depth axis). 



3.54 /xm (Allamandola et al. 1992 ). This feature is ascribed 
to the C-H stretching mode of solid CH3OH. A direct 
comparison with the high mass protostar NGC 7538 : IRS9 
shows that, although the 3.47 fim bands have similar 
shapes, the 3.54 fim feature is absent in Elias 29 (Fig. 
We determine a 3a upper limit to the peak optical depth 
of r(3.54 ^m)<0.036. Scaling with the observed depth and 
column density in NGC 7538 : IRS9 (Brooke et al. |1999D , 
then results in an upper limit to the CH3OH ice column 
density A^(CH30H ice)< 1.5 x 10^^ cm'^, or less than 5% 
of H2O ice toward Ehas 29 (Table |). The other modes 
of CH3OH ice are either much weaker, or are severely 
blended with the strong H2O and silicate bands (e.g. 



the C-0 stret ching mode at 9.7 fj,m; Schutte et al. 1991 



Skinner et al. 1992 ) and thus do not provide better con- 
straints on the CH3OH ice column density. Toward other 
low mass objects, and quiescent dark clouds, low upper 
limits have been set to the CH3OH ice abundance as well. 
The CH3OH ice abundance found in massive protostars is 
generally of the same magnitude (Chiar et al. 1996 ), bu t 
in a few objects significantly larger (Dartois et al. 1999| ), 
than these upper limits. 

3.2.6. The 6.85 /im band 

Elias 29 is the first low mass protostar in which the 
6.85 /im absorption band is detected (Fig. ||). After sub- 
traction of the H2O ice band and the gas phase H2O lines 



(Fig. g), we find that it has a peak optical depth of t ~0.07 
and an integrated optical depth Tint = 7.8 ± 1.6 cm~^. 
When scaled to the H2O ice column density, the strength 
of the 6.85 fj,m band toward Elias 29 is similar to high mass 
protostars (Keane et al. ^000 ). The band profile, e.g. the 
sharp edge at 6.60 ^m, agrees very well with several high 
mass objects, in particular those tracing 'cold' gas and 
dust (NGC 7538 : IRS9, W 33A, GL 989). It clearly devi- 
ates from warmer lines of sight (e.g. S 140 : IRSl; Fig. ||). 
Thus, in this picture, we find that the material responsible 
for the 6.85 /xm band toward Elias 29 is not significantly 
thermally processed. Given the low upper limits to the 
CH3OH ice column density toward Elias 29, only a frac- 
tion of the band, as for high mass objects, can be explained 
by th e C-H bending mode of CH3OH ices (Schutte et al. 
1996). For a detailed band profile analysis and a discus- 



sion o n the origin of the 6.85 /im band, we refer to Keane 
et al. (|2000| ). 



3.2.7. Upper limits to solid CH4, 
OCS, and 'XCN' 



NH3, H2CO, HCOOH, 



Several solid state species have been detected toward lu- 
minous protostars, but are absent toward Elias 29. The 
deformation mode of solid CH4 was detected toward pro- 
tostars, with a peak position at 1303 cm^^ (7.67 /im), 
and a width FWHM=11 cm^^ (Boogert et al. |l996t Dar- 



tois et al. 1998b). For Elias 29 we can exclude this band 
to a peak optical depth of r <0.03, corresponding to 
7V(CH4)/7V(H20)<1.5%. This 3(t upper limit is compa- 
rable to the detection in NGC 7538 : IRS9 (Boogert et al. 
119961 ). 

Solid NH3 was detected by its 9.10 /im inversion mode 
toward NGC 7538 : IRS9 (Lacy et al. [19981). Usin g the 



band strength determined in Kerkhof et al. ( 1999 ), the 
NH3 column density is 13% of H2O ice. Recent detec- 
tions in other highl y obs cured lines of sight give similar 
(W 33A; Gibb et al. |2000|) , or a factor 2 larg er NH3 abun- 
dances (Galactic Center; Chiar et al. 2000|) . To find this 
band in the deep silicate feature of Elias 29 we take the 
same approach as Lacy et al., by fitting a local straight 
line continuum to the wavelength regions 8.52-8.69 and 
9.20-9.55 /im. As a check, we perform the same procedure 
to the ISO-SWS spectrum of NGC 7538 : IRS9 (Fig. §. 
We confirm the detection of Lacy et al., although the 
peak optical depth r ~ 0.16 is a factor 2 lower in our 
case. We ascribe this difference to the calibration uncer- 



tainties of ISO-SWS at this wavelength (Leech |200C| ) . For 
Elias 29, a feature with r ~ 0.06 might be present. How- 
ever, due to the poorly defined long wavelength side of the 
continuum (Fig. |^) and the ISO-SWS calibration uncer- 
tainties, we will assume a conservative upper limit to this 
band of r <0.1. This corresponds to a column density of 
7V(NH3)< 3.5 X 10"cm-2, i.e. A/-(NH3)/iV(H20)<13%. 
Other vibrational bands of NH3 do not provide better 
constraints. The equally strong N-H stretching mode at 
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9.5 8.5 
A (fim) 

Fig. 9. a ISO-SWS spectra of the silicate band region of 
Elias 29 and NGC 7538 : IRS9. The dashed hne is the 
local continuum for the solid NH3 invers ion mode, sim- 
ilar to that defined in Lacy et al. ( 1998 ). b the residu- 
als after continuum subtraction for Elias 29 (solid) and 
NGC 7538 : IRS9 (dotted). 



2.90 i-iia (d'Hendecourt & AUamandola 1986| ) is hidden in 
the steep wing of the 3.0 fim ice band, and there is no sig- 
nificant difference in this region between the laboratory 
spectra of pure H2O ice and the mixture H2O:NH3=10:l 
(Fig. H). A similar problem exists for the N-H deformation 
mode at 6.16 fim, which is hidden in the long w avelen gth 
wing of the H2O bending mode (Keane et al. 2000|) . A 
feature with a peak optical depth of t <0.0 25 would be 
expected here (Sandford & AUamandola 199S ) . When sub- 
tracting water ice and vapor absorption, a weak band with 
an optical depth r — 0.03 perhaps remains present at the 
expected wavelength (Fig. |8|) . Given the other positive and 
negative structure in the spectrum, we regard this also as 
an upper limit, however. 

The H2 0-subtracted 5.0-6.5 fim wavelength region 
(Fig. |8|) does not show the feat ures d e tected toward high 



mass protostars (Schutte et al. 1996 ; 1998 Keane et al. 
2000| ). At 6.25 fj,m (not to confuse with the feature of NH3 
ice at slightly shorter wavelength; see above), a feature 
has been associated with absorption by carbonaceous dust 
(PAH). At 5.83 ^m a broad feature has been assigned to 
the C=0 stretching mode of solid HCOOH, and a narrow 
feature of solid H2CO (Keane et al. 200C). Scaling the fea- 
tures observed toward NGC 7538 : IRS9 to the lower H2O 
ice band column density toward Elias 29, one would ex- 
pect peak optical depths rs.ga = 0.06 and Te,25 = 0.03. Our 
spectra indicate upper limits to these features of t <0.03 
(Fig. ^). Thus, in particular the 5.83 fiui feature toward 
Elias 29 is significantly less pronounced compared to high 
mass protostars. Using the band strengths and typical 
widths given in Keane et al. ( |200C| ), we derive 3cr col- 



umn density upper limits of Af(H2C0) < 6 x 10 cm^ , 
and iV(HCOOH) < 3 x 10^^ cm^^. With abundance upper 
limits of 1-2% with respect to H2O, these aldehydes are 
thus minor ice components. For comparison, toward high 
mass objects it is typically 3% or higher. 

An absorption feature has been detected at 
2042±4 cm~^ (4.90 /im) in lines of sight toward several 
massive protostars (Palumbo et al. 1997). With a width 
FWHM=23±6 cm~^, it has been ascribed to absorption 
by solid DCS. For Elias 29 this feature is not detected 
with a peak optical depth r <0.01 (3(t), corresponding 
to iV(OCS)< 1.5 X 10^5 cm"2 q^. ^ o.05% of H2O ice. 
This upper limit is of the same order of magnitude as the 
detections in W 33A and Mon R2 : IRS2 (Palumbo et al. 



1997) 



Finally, toward several high and low mass protostars a 
feature has been detected at ~2166 cm ~^ (4 .62 ^m) with 
a width FWHM~20 cm"! (Lacy et al. |1984| ; Tegler et al. 



1995 ). This feature is absent in Elias 29, with a peak op- 
tical depth r <0.01 {3a). If this feature is caused by the 
C=N stretching mode in 'XCN', this corresponds to a col- 
umn density iV(XCN)< 6.7 x lO^^ cm'^, or less than 0.2% 
of H2O ice (applying A — 3x 10""'^^ cm molecule""'^; Tegler 
et al. 1995). This is considerably less than the detections 
made toward high mass objects (e.g. W 33A) and several 
low mass objects (Elias 18; L 1551 : IRS5; Tegler et al. 



1995). This feature has not been detected in the quiescent 



regions of the Taurus molecular cloud (Elias 16; Table 
For a more elaborate discussion on this feature, and the 
proposed carriers, we refer to Pendleton et al. ( |1999D . 

3.2.8. Silicates 

The absorption bands of the Si-0 stretching and bending 
modes of silicate dust are prominently present at 9.7 fim 
and 18 /im (Fig. ^). We derive a peak absorption opti- 
cal depth of the 9.7 /im band rg.7 = 1.38 (Fig. which 
is in excellent agreement with the ground-based study of 
Hanner et al. ( 1995| ). It is likely that this is a lower limit, 
since the absorption bands have been partly filled in with 
silicate emission from hot dust near the protostar. Model- 
ing of the 9.7 /im silicate band toward Elias 29, including 
emission and absorption, shows that Tg.7 ranges between 
1.51 and 3.38 for optically thick and thin emission respec- 
tively (Hanner et al. 1995). A better fit is obtained for 
optically thick emission. In contrast, for luminous proto- 
stars optically thin emission has been generally assumed. 
Using the relation T9.7 = 1.4 T9.7(obs) -I- 1.6 (Gillet et al. 
19751 ; Willner et al. |l982|) , yields Tg,r = 3.53 for Ehas 29. 

For these values of r9.7, the visual extinction Ay 
ranges between 28 and 65, assuming the standard rela- 
tion Ay /tq. 7=18. 5 (Roche & Aitken p"984 ). However, these 
limits are likely overestimated (30-50%), because of the 
anomalous extinction curve due to larger grains in the 
p O ph m olecular cloud (Bohlin et al. |1978| ; Martin & Whit- 
tet 1990). Independent extinction determinations, such as 
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Fig. 10. High resolution (i?=2000) 
spectrum of Elias 29 showing many 
gas phase CO lines, and the CO ice 
band at 4.67 fim. The narrow emis- 
sion line at 4.653 /im is Pf /3 of Hi. 
Some of the narrow structure seen 
on the long wavelength side is due 
to lines of hot H2O vapor (see text). 
For comparison we show a low res- 
olution (i?=400) observation, shifted 
along the flux scale for clarity, where 
the gas lines are smeared out over the 
continuum. 



Ay <48 from the H-K broad band color and Ay 
-lis 



<80 



from C O observations (Wilking & Lada 1983 ), do not 
help to solve this issue. Millimeter continuum observa- 
tions (Andre & Montmerlc 1994), and the near- infrared 
J-H color (Greene, priv. comm.), suggest a relatively low 
Ay <30. 

The total hydrogen column density A^h = 
A^(H i)-|-2A^(H2) is closely related to T9.7, and, in 
contrast to the derivation of Ay, the derived A'^h is 
not strongly affected by the large grain size in p Oph. 
Applying standard conversion factors for the diffuse 
ISM (Bohlin et al. |l978t Roche & Aitken |1984[), we 



find A^H = 0.5 — 1.2 x 10 cm^^, depending on the 
applied rg.7. To be consistent with studies of high mass 
protostars, we will assume in the abundance calculations, 
the value corresponding to optically thin silicate emission, 
i.e. the high limit A^h = 1-2 x 10^^ cm"^ (Table |). 



3.3. Gas phase absorption lines 



vations {R — 1500 — 2000). Thus, three parameters are 
varied to fit the observed absorption lines: the column 
density A'', the Doppler parameter &d, and the excitation 
temperature Tex- Reliable column densities can only be de- 
rived if 6d is a priori known, which in many studies (like 
ours) is not the case, since the lines are unresolved. At 
low values of b^, the lines become easily optically thick, 
and much larger column densities are needed to fit the 
observed lines, compared to models with high 6d values, 
and optically thin lines. 

We emphasize that our assumptions of collisional ex- 
citation, and LTE at a single T^x need not be valid. There 
is likely a temperature gradient along the line of sight, 
as expected for a protostellar envelope. The LTE assump- 
tion may not apply for the high rotational levels, which 
have high critical densities. Also, the energy levels may be 
pumped by infrared photons, rather than being collision- 
ally excited. Bearing these caveats in mind, we will here 
focus on deriving CO and H2O gas column densities and 
temperatures using the LTE models. 



The high resolution 4.00-8.50 /im spectrum of Elias 29 
shows an impressive number of narrow absorption lines of 
gaseous CO and H2O (Figs. | and |l|). We determined 
local continuum points by hand and connected these, us- 
ing a smooth cubic spline interpolation. Then the data 
were converted to optical depth scale, and the absorp- 
tion lines were modeled, using the ro- vibrational spectra of 
gaseous CO and H2O described in Helmich ( |1996| ). These 
models assume the gas is in Local Thermodynamic Equi- 
librium (LTE), and has a single excitation temperature 
Tex- The absorption lines have a Voigt profile, and are 
Doppler broadened to a width bu (=FWHM/2\/h^). The 
line oscillator strengths are calculated from the HITRAN 



database (Rothman et al. 1992). Finally, the spectrum is 
convolved with a Gaussian to the resolution of our obscr- 



3.3.1. CO gas 

The 4.4-5.0 /im region shows absorption lines of gas phase 
^^CO, up to rotational quantum number Jiow=33 in the R- 



branch, and Jiow=36 in the P-branch (Fig. 10). The P(l), 
P(2) and R(0) lines are blended with the CO ice band at 
4.67 /im and the H I Pf /3 emission line at 4.653 /xm. For 
all other absorption lines we determined equivalent widths 
to construct a rotation diagram. A rotation diagram gives 
a first impression of the temperature components present 
along the line of sight, as well as their column densities (or 
lower limits for optically thick lines). For technical details 
on constructing such a diagram we refer to Mitchell et al. 
( |1990D , and Boogert et al. ( |1998D . The equivalent widths 



were converted to column densities, using the oscillator 
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Fig. 11. a rotation diagram of the ^^CO lines detected 
toward Elias 29 showing the presence of hot and cold gas 
along the line of sight by the different slopes at high and 
low rotational levels, b xt contour diagram of model fits 
to the observed ro-vibrational spectrum of the R-branch 
of gaseous CO toward Elias 29. xt values are shown for 
the temperature Tex versus CO column density N at con- 
stant velocity broadenings 6d of 5 km s~^ (dotted) and 
10 km s~-^ (solid). We only show models that provide ac- 
ceptable fits to the data, i.e. xl < 4- 



strengths of Goorvitch (|994|). For ^^CO (Fig. |TT]), we find 
two regimes with very different slopes, corresponding to 
temperatures T^t = 90 ± 45 K and T^t = 1100 ± 300 K 
respectively (with 3a errors). However, the slopes of the 
R- and P-branch lines of the hot component are different 
(Fig. |ll|), resulting in T^t = f700 ± 420 K when fitting 
to the P-branch lines only. A possible explanation is that 
CO is excited by continuum photons rather than collisions. 
The rising continuum may lead to a higher T^ot for the P- 
branch with respect to the R-branch. This effect becomes 



Fig. 12. a Rotation diagram of the ^^CO lines detected to- 
ward Elias 29. Circles are P-branch lines, and triangles are 
R-branch lines. Open symbols refer to ^'^CO lines heavily 
blended with ^^CO lines, and are not used to determine 
the physical parameters. The straight line indicates the 
best fit, with a gas temperature Trot = 85±57 K, and a col- 
umn density N(^^ CO) ^{1.1 ±0.2) x 10^"^ cm'^. b xl con- 
tour diagram of model fits to the observed ro-vibrational 
spectrum of gaseous '^'^CO toward Elias 29. xt values are 
shown for the temperature Tox versus CO column density 
N for constant velocity broadenings 6d=2.5 km s^^ and 
6d=10 km s~^. Only acceptable fits to the data, having 
xt < 3.5, are shown. 



stronger when the photons released after de-excitation of 
R-branch levels are re-absorbed in P-branch levels. Ra- 
diative excitation has also been used to explain the H2O 
ro-vibrational spectrum toward Orion BN/KL (Gonzalez- 
Alfonso et al. |1998D. The fact that the H2O P-branch hues 



are seen in emission for Orion BN/KL, rather than in ab- 
sorption as for CO (and H2O; Sect. 3.3.2) toward Elias 29, 



12 



A.C.A. Boogert et al.: The low mass protostar Elias 29 



may reflect a different density gradient toward Elias 29, 
such that the photons are not able to escape the enve- 
lope. Additionally, coUisional excitation in shocks may be 
of less importance in Elias 29 compared to Orion BN/KL. 
A more careful analysis is needed to discriminate between 
the radiative and coUisional excitation models, and alter- 
native explanations, such as non-LTE effects. 

The column densities that we derive from the ab- 
scissa in the rotation diagram are A^(CO)= 1.7 x 10^^ and 
7V(C0)= 3.5 X lO^'^ cm-2 for the cold and hot CO compo- 
nents respectively. To better constrain the column densi- 
ties and derive more reliable temperatures, one has to take 
into account optical depth effects, using the LTE model 
spectra discussed above. We chose to fit to the frequency 
range 2170-2290 cm~^ (Jiow > 7 in R- branch), thus mini- 
mizing the contribution from the cold CO component and 
contamination by ^^CO lines (see below). We find that 
good fits to these high R-branch lines are obtained only 
for line widths 6d> 3 km s^^. Sub- millimeter emission line 
studies indicate 6d=3.6 km s^^ for CO J= 6^5, but 
much lower values of 6d=1.2 km s~^ for C^*0 J= 1-^0 
and CS J= 5 — *■ 4 (Boogert, Hogerheijde, et al., in 
prep.). Indeed, studies of other sources have shown that, 
as a rule, infrared absorption lines are broader than sub- 



Table 2. Gas phase ^^CO and H2O column densities, de- 
rived with various methods 



millimeter emission lines (van der Tak et al. 1999). Fig- 
ure |l^ shows the xt contour diagram of temperature ver- 
sus column density for two values of the line width 5d=5, 
and 6d=10 km s^^. The best fitting models have temper- 
atures Tex = 1100±400 K, in good agreement with the 
rotation diagram. At &d=10 km the column density is 
well constrained to A^(C0)=(1.3± 0.5)xl0is cm'^, which 
is a factor of 3 larger compared to that derived from the 
rotation diagram. Thus at 6d = 10 km s"'^ the lines are 
still somewhat optically thick. At lower 6d=5 km s~^, the 
lines become very optically thick, and the column density 
is poorly constrained. Although the best fits with x^< 3 
have iV(CO)=(8±4) x lO^^ cm'^ at Tex = 650±150 K, rea- 
sonable fits are obtained at any N{CO)> 2 x 10^* cm^^ 
for this hot CO gas. 

Several ^^CO lines can be seen in between the ^^CO 
P-branch lines (Fig. |l^). At the resolution of our observa- 
tions, the blending with the ^^CO lines hinders analyzing 
the much weaker ^'^CO lines. But using several well sepa- 
rated lines, we were able to construct a rotation diagram 
(Fig. |l2|). We find that they result from cold gas at Tj-ot = 
85±57 K (3(7 error), in good agreement with the cold ^^CO 
gas temperature. In the optical thin case, the column den- 
sity of this cold component is N(^^ CO) ={1.1 ±0.2) x 10^^ 
cm~^. However, the detected ^"^CO lines could still be op- 
tically thick. Therefore, we also modeled the ^'^CO spec- 
trum, and determine the xt after subtraction of a good 
fitting hot i^CO gas model (Figs. ^ and O). In the op- 
tically thick case, such as for 6d=2.5 km s , the column 
density can have a wide range N(^^C0)={2 ± 1.3) x 10^'' 
cm~^. Using the isotope abundance ratio ""^^00/^^00=80 
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sity is thus N(^^ CO) ={16 ±10) x lO^^ cm'^. There is also 
evidence for ^"^00 lines of warm gas (Jiow > 9), but at low 
significance (< 2(j) and no reliable temperature or column 
density could be derived. 

We conclude that the CO gas along the line of sight 
consists of two temperature components, Tj-ot = 90±45 K 
and Trot = 1100±300 K. The column density of both com- 
ponents depends highly on the assumed line optical thick- 
ness (Table ||). Until the intrinsic line width is directly 
observed by very high spectral resolution observations, 
we can only give a lower limit of A^(CO-hot)> 2 x lO""^* 
cm~^, while iV(CO-cold) is not well constrained, i.e. 



(16 ±10) x 10^8 cm 



-2. Given that Nh = 1.2 x lO^^ cm^^ 
toward Elias 29, a total gas phase CO column density 
Af(C0)=12 X 10^^ cm~^ is expected, assuming that most 
of the gas along the line of sight is molecular and the con- 
version factor iV(H2)/N(CO)=5000 applies (Lacy et al. 
Then, the ratio of hot to cold CO gas along the 



1994) 



line of sight must be at least 0.2. 



3.3.2. H2O gas 



(Boogert et al. 200C), the inferred cold ^^CO column dcn- 



We compare the numerous narrow absorption lines de- 
tected in the 5-7.3 fj,m spectral region of Elias 29 with 
model spectra of H2O vapor at various physical conditions 
(Fig. H). Clearly, the many lines observed at wavelengths 
longer than ~6.55 /xm are explained by II2O vapor at a 
high temperature (Tex >100 K). On the other hand, the 
relative weakness of the lines observed in the range 6.55- 
6.65 /im imposes a strict upper limit to the temperature 
of this hot gas (Tex <1000 K). To further constrain the 
gas temperature, and the II2O vapor column density, we 
determined the xt for a large number of models. Rea- 
sonable fits to the full 5-7.3 /im range are obtained for 
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Fig. 13. The P branch of gas phase CO observed in Ehas 29 (top) compared with a weU fitting ^^CO gas model at 
Tex = 750 K (A'' = 5 X 10^^ cm~^, 6 = 5 km s~^). The bottom panel shows the residual after subtraction of the -^^CO 
gas model, which contains lines of cold ^^CO gas. For comparison, a ^^CO model is plotted (Tgx = 50 K, A'' = 2 x 10^'' 
cm~^, b = 2.5 km s~^). 
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Fig. 14. Optical depth plot of gaseous H2O lines observed toward Elias 29 (top), compared with model spectra at 
various temperatures. Column densities are iV = 5 x 10^^ cm~^ for all models, and 6d=5.0 km s~^ (Tex = 50 K models 
have 6d=2.5 km s~^). The single component T = 300 K and the two component T = 50 + 500 K models give equally 
good fits. The models with just cold (T = 50 K) or very hot (T = 1000 K) gas clearly do not fit the data. 



temperatures of Tex = 350 ± 200 K. The column density is depths (&d>5). For narrower lines the column density can 
constrained to A''=(7±4) x 10^^ cm~^ for low line optical be an order of magnitude larger. 
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Table 3. Line-of-sight averaged solid and gas phase abundances {N/N-r in units of 10^^) 



Species 


Dense Cloud" 


N7538/9 


W 33A 


Elias 29 


GL 2591 


Refs.'' 


H20-ice 


64 


50 


39-143 


28 (8) 


10 


[1], [2], [3,4], [5], [6] 


-gas'^ 


<1 


< 3.1 


<3.6 


>3 


24 (3) 


[7], [8], [8], [5], [8] 


co-ice'' [total] 


17(1) 


8(1) 


3.2 (1.8) 


1.4 (0.2) 


<< 1 


[9], [9], [9], [5], [10] 


-ice [apolar] 


14 (1) 


7 (0.5) 


0.8 (0.1) 


1.2 (0.2) 


<< 1 


[9], [9], [9], [5], [10] 


-gas 


< 10 


91 (51) 


143 (32) 


> 67'' 


113 (15) 


[11], [12], [12], [5], [12] 


C02-ice 


12 (3) 


10 (1) 


5.2 (0.5) 


5.4 (0.5) 


0.9 (0.1) 


[13], [13], [13], [13], [13] 

L J ! L J ' L J ' L J ' L J 


-gas 




0.05 (0.01) 


0.08 (0.02) 


<0.06 


0.15 (0.03) 


[8], [8], [8], [8] 


NHs-ice 




7.6 


6.1 


< 3.0 




[14], [3], [5] 


CHgOH-ice 


< 1.8 


2.0 


7 


< 1.3 


4(2) 


[15], [15], [15], [5], [16] 


HaCO-ice 




1.9 


2.5 


< 0.5 




[4], [4], [5] 


HCOOH-ice 




1.1 


0.6 


< 0.3 




[4], [4], [5] 


CH4-ice 




0.8 


0.6 


< 0.4 




[17], [17], [5] 


OCS-ice 


< 0.13 




0.07 


< 0.02 




[18], [18], [5] 


'XCN'-ice 


< 1.3 


1 (0.3) 


3.6 (1) 


< 0.06 




[19], [19], [19], [5] 


A'"h [10^^ cm-^f 


0.39 


1.6 


2.8 


1.2 


1.7 


[20], [20], [20], [5], [21] 



Ice abundances are for Taurus dense cloud toward Elias 16. Gas p hase abundances are fo r p Oph cloud 



References from left to ri ght for each column: [1] Chiar et al , 
4] Keane et al. |20 0C|; [5 ] this work; [6] Smith et al 



200C 



al. 2000; [9] Chiar et al 



1999 



[13] Gerakines et al 
et al. |1998|; [18] Palumbo et al 



1998 



10] van Di shoec k et al. 1996; [11] Caux et al 
[15 



[14] Lacy e t al 



119971; 

T 



1998 



1989 



1995 



[7] Liseau & Olofs son 



[19] Tegler et al 



Chiar et al, 
[20' 



1995 



1996 



[2] Schutte et a l. [l996| ; [3] Gibb et al. 

[8] Boonman et 



199£ 



1999; [12] Mitche ll et al 



[16] Schutte et al 
Tielens et al 
for cold gas 



pjgi 



1990 



^ ^, 17] Boogert 

199l| ; [21] this work. Fig. | 



All models assume 6d=5 km s^ for hot gas and 6d=2.5 km s' 
Total CO ice abundance given as well as the abundance present in the polar and apolar ice components 
along the line of sight 

" Determined from the 9.7 pm silicate band 



In a second approach, we test whether both hot and 
cold H2O vapor components could be present along the 
line of sight, much like the hot and cold CO compo- 
nents. We fitted the regions 5.5-5.8 and 6.55-7.3 /xm, 
which do not contain lines from the lowest rotational 
levels and thus are particularly sensitive to warm H2O 
vapor along the line of sight (Helmich et al. 1996; Dar- 
tois et al. 1998b). The excitation temperature of this gas 
is Tox — 500 ± 300 K, with column densities similar to 
that of the single component model. In the high temper- 
ature regime (Tex >500 K), the modeled line depths in 
the 6.0-6.5 fim region, tracing colder gas, are significantly 
underestimated. To determine the temperature and col- 
umn density of this possible cold component, we fitted 
the sum of a good fitting hot gas model (Tex =500 K, 
N = 5x 10^^ cm~^, 6d=5.0 km s~^) and a grid of models 
at a wide range of physical conditions to the spectrum of 
Elias 29. Thus, here we assume that the lines of the hot 
and cold gas have different radial velocities and the optical 
depth spectra can simply be added. We find that indeed 
a significant amount of 'cold' II2O vapor, at Tox < 200 K 
may be present (Fig. |l|). At Tex < 100 K the column 
density exceeds the assumed hot H2O column density of 

= 5 X 10^^ cm^^. For a line width of 60=5. km s^^, we 



find that N <lx 10^^ cm"^. At 6d=2.5 km s"\ the col- 
umn density of this cold II2O gas cannot be constrained. 

To summarize, the lines in the 5-7.3 /im range are rea- 
sonably fitted with II2O models at T^x ~ 350 ±200 K, and 
A^=(7±4) X 10^^ cm^^ at low fine optical depths. For nar- 
rower lines (&]3<5 km s^^), the column density can be an 
order of magnitude larger. In accordance with the gaseous 
CO along the line of sight, equally good fits are obtained 
with a two component model, where the cool component 
has Tox < 200 K, and the warmer component Tox >500 K. 
The cool component is then at least as abundant as the 
warm H2O gas. 



4. Discussion 



.1. Gas and solid state abundances 



We have calculated line of sight averaged gas and solid 
state abundances toward Elias 29, by dividing the col- 
umn densities derived in this paper over the total hydro- 
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gen column density Nh — 1-2 x 10 cm^ (Sect. ^ 
We compare these abundances with a sample of sight- 
lines, spanning the range from dark cloud core to fairly 
evolved protostars (Table |3|) . As a tracer of ices in quies- 
cent dark cloud material, we chose the object Elias 16, an 
evolved star by chance located behind the Taurus molec- 



ular cloud (e.g. Whittet et al. 1998). Gas phase CO and 
H2O abundances in dense clouds were taken from ISO- 
LWS studies (Caux et al. |1999| ; Liseau & Olofsson pJ99| ) . 
The least evolved protostar in our comparison sample is 
NGC 7538 : IRS9. The infrared spectrum of this deeply 
embedded object is characterized by cold ice (Whittet et 
al. 1996| ), and the gas phase temperatures and abundances 
indicate a very mode st hot core (Mitchell et al. |l990| ; 
Boonman et al. ^OOOl) . W 33A is more embedded than 
NGC 7538 : IRS9, but does have a significant amo unt of 
warm gas along the li ne of sight (Mitchell e t al 
Lahuis & van Dishoeck 2000; Boonman et al 



199C 



200c ), an d 



has a lower abundance of volatile ices (Tielens et al. 1991 ) 
The most evolved object in our sample is GL 2591. It is 
a typical high mass hot core source, with low ice abun- 
dances and high gas temperatures. All these protostars 
are associated with infrared reflection nebulae, and have 
well developed high velocity molecular outflows (Mitchell 
et al. |1991| ; Bontemps et al. |l996t van der Tak et al. |2OOO0 . 
Finally, it is important to note that all the comparison 
protostars are at least three orders of magnitude more lu- 
minous than Elias 29. This allows an investigation of the 
effect of low and high mass star formation on the molecular 
envelopes. An extensive comparison with low luminosity 
embedded objects is at present not possible, because their 
infrared gas and solid state characteristics have not been 
studied in such great detail. 

The H2O and CO ice abundances decrease for the se- 
quence of quiescent dense cloud to NGC 7538 : IRS9, 
W 33A and GL 2591 (Table |). At the same time, the gas 
phase H2O abundance, the gas phase CO and H2O tem- 
peratures, as well as the gas-to-solid ratios (Table in- 
crease for these objects. All these effects can be explained 
by evaporation of the ice mantles and heating of the hot 
core. It has been suggested that the observed H2O gas 
may also have been newly formed by reactions of atomic 
O and II2 in warm conditions (T > 200 K) in the cen- 
tral hot core or in s hocks created by the outflow (e.g., 
van Dishoeck & Blake 1998 ). However, the total (gas plus 
ice) H2O abundance decreases for the more evolved ob- 
jects, indicating that II2O is destro yed rather than being 
newly formed (van Dishoeck |199^ ). The low gas phase 
CO2 abundance in all sources indicates that this molecule 
is destroyed even more effici ently after evaporation from 
the g rains (Boonman et al. 200C; Charnley & Kaufman 
|2000|). 



In the proposed heating sequence, Elias 29 is placed af- 
ter W 33A, and before GL 2591. However, the various ice 
band profiles (H2O, CO, CO2, and 6.85 /im) in Ehas 29, 
indicate little thermal processing, resembling very much 
NGC 7538 : IRS9, rather than W 33A or GL 2591. The 
combination of high gas phase abundances and tempera- 
tures, together with a lack of signatures of thermal pro- 
cessing in the ice bands, as seen in Elias 29, is remarkable 
and is not seen in high mass protostars. Geometric effects 
may play an important role in the evolution of molecular 



envelopes around low mass protostars (see Sect. 4.2). 

Whereas thermal evaporation can explain the abun- 
dance variations of volatiles such as H2O, CO, CO2, NH3, 
and CH4, other mechanisms are needed to explain the vari- 
ations of solid CH3OH, and XCN abundances among the 
sources in our sample (Table ||). It has been widely consid- 
ered that XCN molecules are formed by energetic process- 
ing of icy grain mantles by stellar or cosmic ray induced 
far-ultraviolet radiation, or by bombardment with highly 
energ etic p articles (e.g. Lacy et al. 1984 , Grim & Green- 
berg |l987t AUamandola et al. |l988[ ). The high CH3OH 
abundances toward sources with deep XCN bands, and 
the apparent absence of CH3OH toward low mass proto- 
stars and dark clouds might suggest that the energetics of 
nearb y ma ssive stars is needed to produce CH3OH (Gibb 
et al. |200C|). 



Table 4. Gas-to-solid state column density ratios 



Object 


CO" 


H2O" 


CO2" 


r-p b 

J- warm 


Dark CI. 


<1 


<0.02 






N7538/9 


12 (6) 


<0.05 


0.005 


180 (40) 


W 33A 


45 (22) 


<0.11 


0.015 


120 (20) 


Elias 29 


>53 


>0.23 


<0.011 


1000 (500) 


GL 2591 


>400 


2.4 


0.17 


1000 (200) 



" Determined from Table ^ 
Temperature of warm CO gas (Mitchell et al. |l990[) 



^ For actual, local, ab unda nces in high mass protostars we 
refer to Boonman et al. 200C 



4-.2. The structure of Elias 29 

The variety of dust, gas and ice absorption and emis- 
sion components presented here, and in the literature, al- 
lows us to construct an overall view of the structure of 
Elias 29. The scale on which the detected hot CO gas is 
present can be constrained when one assumes that the 
pure rotational high-J CO emission lines detected toward 
Elias 29 with ISO-LWS (Ceccarelli et al., in prep.) are 
emitted by the same hot gas. We fit these observed line 
fluxes, by assuming spontaneous, optically thin emission 
from an LTE level distribution, and leaving the size of 
the emitting region as a free parameter. For the range of 
column densities and temperatures found to fit the CO ab- 
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sorption lines (Fig. , we find diameters in the range 85- 
225 AU. Thus, the observed hot CO gas may be present in 
a hot core region with the size of a circumstellar disk. The 
gas could be concentrated in a high density photospheric 
layer above the disk. To sufficiently heat it by radiation 
from the central star, the disk needs to flare outwards, 
rather than being flat (e.g. Chiang & Goldreich 1997). We 



cannot exclude however that the gas is present more uni- 
formly in the hot core, at lower densities. It might then 
also be partly heated by shocks from the outflow close to 
the star. A more detailed modeling of the CO emission 
lines, including departures from LTE, optical depth cor- 
rections, and taking into account excitation by radiative 
pumping, are needed to further confine the location of the 
gas phase CO and H2O components (Ceccarelli et al., in 
prep.). 

The lack of signatures of thermal processing in the 
ice bands, locates the ice in a region shielded from the 
central heating source. The ices could be present in a 
foreground cloud, an extended envelope, or a circumstel- 
lar disk seen close to edge-on. Millimeter wave continuum 
observations indicate an extended envelope, concentrated 
on the infr ared s ource (FWHM=17"; 2600 AU; Andre & 

199^ ). For a wide range of 



Montmerle 1994; Motte et al 



power law model fits to the far-infrared SED, Andre & 
Montmerle find that the envelope mass is 0.1 Mq, with a 
volume averaged dust temperature of typically T =35 K. 
This temperature is high enough to evaporate the most 
volatile, apolar ices, but too low to induce ice crystalliza- 
tion. Hence, indeed the observed H2O, CO2, and probably 
"6.85" ^m ices could be associated with this extended en- 
velope. Some of the apolar CO ice has evaporated in the 
envelope after the formation of the low mass protostar, 
as indicated by the significantly lower CO/H2O ice ratio 
toward Elias 29, compared to other sight-lines with little 
thermal processing in the ices (NGC 7538 : IRS9, Elias 16; 
Table ^). In this picture, the detected apolar CO ice could 
thus be spatially separate from the other ices, perhaps in 
foreground clouds, or well shielded in a very cold disk. 

Knowledge of the source structure is essential to in- 
terpret the observed solid and gas phase species. For ex- 
ample, if the ice is present in the disk, rather than in the 
envelope, we must see the disk in a near edge-on config- 
uration. Is there independent evidence for the presence 
of a disk surrounding Elias 29 and what would be its 
orientation? The most direct view is provided by lunar 
occultation observations. A central object with diameter 
of ~1 AU emits 90% of the 2.2 /im continuum emission 
(Simon et al. 1987| ). The remaining 10% comes primarily 
from an object of 60 AU in diameter, which could be the 
hot part of a disk (T ~ 1000 K). The strongest spectro- 
scopic disk indicator would be the presence of emission 
or absorption of vibrational overtone band heads of CO 
(e.g. Carr |1989| ; Najita et al. |1996D . The 2.0-2.5 nui spec- 
trum of Elias 29 does not show these features, in contrast 
to other protostars in p Oph, such as WL 16 (Greene & 



Lada 1996 ). However, the absence of CO overtone bands 
does not prove the absence of an (inner) disk (Calvet et 
al. 1991). For example, the Herbig Ac object AB Aur does 
not have detected CO overtone bands, while high spatial 
resolution radio continuum and emission line observations 
provide strong evidence for the presence of a circumstellar 
disk around this object (Mannings & Sargent 1997). 

AB Aur is an interesting comparison source, since it 
has the same luminosity as Elias 29 (^--^40 Lq), and the 
SEDs of both objects are remarkably similar (Fig. H). The 
flatness of the SED in AB Aur is well reproduced in flar- 
ing disk models, where the dust in the outer parts of the 
disk is more efficiently heated than in flat disks (Chiang & 
Goldreich 1997). The disk is optically thick up to 100 ^m, 
and becomes optically thin at longer wavelengths where 



the SED drops steeply (e.g. van den Ancker et al. 2000). 
The similarity of the SEDs does however not necessarily 
imply that Elias 29 is dominated by an optically thick disk 
as well. A flat SED could also be produced by the enve- 
lope, if it has a shallow po wer la w density profile (index 
'--^ 0.5; Andre & Montmerle 1994 ). This density profile is 
remarkabl y fiat compared to h igh ma ss protostars (van der 



Tak et al. 2000; Dartois et al. 1998a ), and ot her low mass 



protostars (e.g. Hogerheijde & Sandell 2000 ). Finally, flat 
energy distributions are also created by the combination of 
a disk and envelope. Here, the he ated e nvelope irradiates 
the outer parts of the disk (Natta [l993 ). 

Without direct high resolution imaging, it is difficult 
to discriminate between these models. Assuming a given 
model however, the present observations put some con- 
straints. In the disk scenario, its orientation would have 
to be closer to edge-on than face-on to explain the ab- 
sorption line spectrum of Elias 29 (Chiang & Goldreich 
1999 ). In these models, an inclination larger than ~70° 
can however be excluded, because this would give an SED 
that peaks in the far-infrared, in contrast to what is ob- 
served for Elias 29. Also, if the disk were edge-on, a higher 
absorbing column, perhaps an order of magnitude larger 
than the observed A^h ~ 1-2 x lO^^ cm'^ (Sect. U) would 



be expected (Sekimoto et al. 1997). An independent mea- 
sure for A^H and the disk inclination is provided by the 
hard X-ray flux and spectrum, arising from hot gas in the 
magnetosphere. For Ehas 29, a high A^h ~ 2 x 10^^ cm~^ is 
observed during X-ray flares, but A'h is a factor of 5 lower 



in quiescent phases (Kamata et al. 1997). Perhaps the X- 
ray flares are formed low in the magnetosphere, and in the 
relatively high inclination of the disk, they trace higher 
column densities compared to X-rays formed in quiescent 
phases higher in the magnetosphere. 

5. Conclusions and future work 



The 1.2-195 /im spectrum of the low mass protostellar 
object Elias 29 in the p Ophiuchi molecular cloud shows a 
wealth of absorption lines of gas and solid state molecules. 
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Hot CO and H2O gas are detected (Tex >300 K) at rather 
high abundances, on scales of not more than a few hun- 
dred AU. The ice abundances are relatively low. In this 
respect, Elias 29 resembles luminous protostars with sig- 
nificantly heated cores, such as GL 2591. However, none of 
the many ice bands that are detected, i.e. from H2O, CO, 
CO2, and the 6.85 /^m band, shows outspoken signs of 
thermal processing. Again in comparison with luminous 
protostars, Elias 29 now resembles less evolved objects, 
such as NGC 7538 : IRS9. Our combined gas and solid 
state analysis thus shows that high and low mass pro- 
tostars heat their molecular envelopes in different ways. 
This may be related to their different structure, such as 
the presence of a circumstellar disk in low mass protostars. 
The hot gas of Elias 29 could be present on the surface 
of a flaring disk, which is efficiently heated by the central 
star. The ices toward Elias 29 must be well shielded in a 
circumstellar disk seen close to edge-on, or far away in the 
envelope. 

Does this imply that in general the ices in the disks 
or outer envelopes of low mass protostars remain unal- 
tered, both in composition and structure, during the pro- 
cess of star formation? Are these ices the building blocks 
of the early solar system and are they preserved in present 
day observed cometary nuclei? To date, no Class I proto- 
star has been found with strong signs of crystalline ices 
(Boogert et al. 2000D . On the other hand, the presence of 
crystalline ices and silicates has been reported in several 
isolated, less embedded Herbig AeBe objects (Malfait et 
al. 1999| ). This research needs to be extended to a larger 
sample of low mass protostars, in a range of evolutionary 
stages and luminosities. Furthermore, it is essential for the 
interpretation of the gas and solid state characteristics to- 
ward Elias 29 that the presence of a circumstellar disk, 
and its inclination are determined by future high spatial 
resolution infrared or millimeter continuum observations. 
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